clear all
%load ParEst_eh1_time1_Mex ParEst_eh1_time1_Mex
load ParEst_eh1_time0_Mex ParEst_eh1_time0_Mex
Chute = [-1.4, ParEst_eh1_time0_Mex(17:20)]; % gamma, alphaf_2, betaf_2, alphai_2, betai_2 
%Chute = [ParEst_eh1_time1_Mex, ParEst_eh1_time0_Mex(17:20)]; % gamma, alphaf_2, betaf_2, alphai_2, betai_2 
HandleLo = @Principal_A2;
LoBo = [-2, 0.1, 0.1, 0.1, 0.1];
HiBo = [2, 10, 10, 10, 10];
                 options = saoptimset( ...
                     'PlotFcn',{@saplotbestf,@saplottemperature,@saplotf,@saplotbestx,@saplotx},...
                     'TolFun',3000,'ReannealInterval',4, ...
                     'InitialTemperature', 10*10^2*[4 10 10 10 10], ...
                     'TemperatureFcn', @temperaturefast, 'StallIterLimit', 1500);
[x,fval,exitflag,output]  = simulannealbnd(HandleLo, Chute, LoBo, HiBo, options)
ParEst_eh1_time1_Mex = x;
save ParEst_eh1_time1_Mex ParEst_eh1_time1_Mex